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ABSTRACT 

We present early high resolution spectroscopic observations of the nova VI369 Cen. We have detected an 
absorption feature at 6695.6 A that we have identified as blue-shifted ’Li I A6708 A. The absorption line, 
moving at -550 km/s, was observed in five high-resolution spectra of the nova obtained at different epochs. On 
the basis of the intensity of this absorption line we infer that a single nova outburst can inject in the Galaxy 
Mu = 0.3 - 4.8 xl0“'° Mq. Given the current estimates of Galactic nova rate, this amount is sufficient to 
explain the puzzling origin of the overabundance of Lithium observed in young star populations. 

Keywords: novae, cataclysmic variables — Galaxy: abundances 


1. INTRODUCTION 

The light elements Deuterium, 3-Helium, 4-Helium, and 
7-Lithium are synthesized in non-negligible amount s during 
the first few minutes of the initial cosmic expansion ( |Kolb & 
|Turner|1990[ ). For Big-Bang Nucleosynthesis in the standard 
model of cosmology and particle physics (SBBN), the pre¬ 
dicted abundances for these elements depends only on one pa¬ 
rameter, the baryon-to-photon density ratio. Recently, Planck 
data have determined the baryon density to excellent preci¬ 
sion, leading to a primordial lithium abundance in the range 
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A{Li) = 2.66-2.7^ ( Coc et ar]|2014| l. This value is signif¬ 
icantly larger than A{Li) ~ 2.1 -2.3 obtained for old metal- 
poor ([Fe/H] < —1.4) halo stars, whose distribution in the 
Lithium abundance - metallicity diagr am is almost flat, and 
defines the so-called ’’Spite plateau” ([Spite & Spite] 1982 


Bonifacio et al.|2007|l. These stars were long thought to share 


a common, primordial, Li abundance. The discrepancy may 
be explained by diffusion and turbulent mixing in the stel¬ 
lar interiors ([Korn et al.||2006|l, which corrodes the primor- 
dial Li abundance , and/or by a non-standard BBN scenario 
(locco et al. 2009|l. Yet, another puzzle exists. The abundance 


‘ A(Li) = logw(Nn^jN„)+ 12 
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of Lithium observed in the upper envelope of young metal- 
rich ([Fe/H]^ > -1.4) stars traces a growth of the Li abun¬ 
dance from me Spit e plateau value to the meteoritic value, 
A{Li) = 3.26 ±0.05 ( [Lodders et al.||200^ and even higher, 
which clearly indicates mat limium enrichment mechanisms 
must occur on Galactic scales. 

In recent decades several potential Li producers have been 
proposed on theoretical grounds. Among the possible lithium 
farms in the Galaxy, the most plausible astrophysical so urces 
are re presente d by asymptotic giant branch (AGB) stars ( |Iben| 
fT973]l, novae (jStarrfield et al.|19781l and Galactic cosmic ray 
spallation (via fragmentation of material due to impact of ac- 


celerated protons, leading to expulsion of nucleons. 

Lemoine 

et al.|l998ll. Galactic chemical evolution models (|] 

)’Antona 

& Matteucci|1991 Romano et al.|20011 Travaglio et al.|2001 

Prantzos|2012|l, show a convincing match with observations. 


when all lithium producers described above are included. 

The main channel for thermonuclear production of lithium 
inside stars is based on a particular sequence of events, known 
as hot bottom burning and beryllium transport mechanism 
(Cameron & FowI^|1971|). At first there is formation of 


Be via the reaction Tie +a, which requires large temper¬ 
atures T > IQ^K, and occurs in the deep interior of stars. 
Outward convection mechanisms concurrently must trans¬ 
port the newly formed ^Be to external cooler regions, where 
T « lO^/f. Here, beryllium is able to capture electrons (free 
and bound), giving rise to ^Li with the additional formation 


of 0.86 MeV neutr inos. More recent computations ( Ventura 
& D’Antona|2010|l suggest that a large abundance ofTTcan 


thus be achieved in AGB stars, particularly with mass larger 
than 7Mq. However, large lithium yields and, consequently, 
a significant enrichment in Li of the surrounding interstellar 
medium, can be obtained from these sources only if exten¬ 
sive mass l oss is associated with the phases of maximum Li 
production ( jRomano et al.|200T| l. 

In nova systems the abundance of lithium is related to the 
thermonuclear runaway (TNR) scenario, which accounts for 
the physical explanation of ^Li and its ejection from the bi¬ 
nary system. In TNRs the expanding gas reaches velocities 
of the order of 400-4000 km/s, which results in escape from 
the system, and enrichment of the interstel lar medium with 
eleme nts synthesized in the nova outburst Pose & Hernanz| 
|1998| l. This scenario remains uncertain because of the lack 
of direct detection of ^Li during nova outburst, althoug h its 
presence was proposed long ago ( Starrfield et al.|[r978 l. In¬ 


deed, ^Li is easily destroyed by proton fusion at temperatures 
greater than 2.6 x 10^ K, giving rise to two helium atoms. For 
this reason, ^Li detection in stars must focus on the external 
regions, where the temperature is sufficiently cool to allow 
non-depletion of ^Li. 

In the literature there have been reported only upper lim¬ 
its to the abundance o f ^Li in selected novae, e.g., HR Del, 
IV Cep, and NQ Vul (Friedjung 1979|l, while direct detec¬ 


tions of lithium in symbiotic novae (e.g. T CrB, RS Oph and 
V407 Cyg) are characterized by velocities of a few dozens 
km/s, which implies that lithium is not associated with the 


WD ejecta, but rather originates in the secondary giant star 

(jWallerstein et al.|2008 

IBrandi et al.|2009l|Shore et al.|201 

Recently Tajitsu et a 
^Be features in the nea 
Delphini 2013). Indee 

1 . ( 
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3, 

2015 1 have announced detection of 
JV late spectra of V339 Del (Nova 
the Cameron-Powler mechanism is 




expected to occur during the TNR (jStarrfield et al.||1978|), 
leading to an over-production of beryllium at the epoch of 
the nova outburst, with a consequent decay to lithium after 
^ 50 days. The amount of Be observed is larger than typi¬ 
cal values deduced from theoretical p redictions for CO novae 
from thermonuclear runaway models (|Truran|19811l represent¬ 
ing the very early phases of the nova outburst. 

2. OBSERVATIONS OE ’LITHIUM I 6708 A 

We report here the detection of ^Li I 6708 A in the early 
(first three weeks) high-resolution spectra of the “slow” (f 2 = 
40 ± 1 days, where f 2 corresponds to the time during which 
the nova declines by 2 magnitudes) nova VI369 Cen (Nova 
Cen 2013), when the nova was in the optically thick phas^ 
Optical high-resolution spectral observations started 4 days 
after the initial outburst. Early spectral data were obtained on 
days 7, 13, and 21 with the Fiber-fed Extended Range Opti¬ 
cal Spectrograph (EEROS - R ^ 48 000) mounted on th e MPG 
2.2m telescope located on La Silla (|Kaufer et al.|1999|l, and on 
days 4, 11, 16, and 18 with the PUC High Echelle Resolution 
Optical Spectrograph (PUCHEROS - R ^ 20000) mounted 
on the ESO 0.5m telescope located at the Pontificia Uni versi- 
dad Catolica Observatory in Santiago ( |Vanzi et al.j2012 l. The 
spectra are all characterized by the presence of bright Balmer 
and typical Ee II emission lines, which su ggests that th e nova 
was engulfed in its ’’iron curtain” phase ( |Shore|[T012| l. The 
presence of many absorption features in the range 3700-4600 
A in the first two weeks (days 4, 7, 11, and 13) complicates 
the identification of the most common transitions detected in 
the optically thick phase of novae, which are identified from 
their respective P-Cygni absorptions. We identify multiple 
expanding velocities for each transition : in the first week 
(days 4 and 7) we measure from the Balmer lines, 017773-7, 
8446 A, and Ee II (multiplet 42) lines two expanding systems 
with mean velocities of ~ -550 and -1350 (with a maximum 
value of -1400) km/s, while for Na I we see only the system 
expanding at Vexp ~ -550 km/s. In the second week (days 11 
and 13) we identify two prominent expanding components at 
Vexp ~ -550 and -1100 km/s from P-Cygni profiles of all main 
transitions (Balmer, O I, Ee II and Na ID), while an additional 
faster component is observed at Vexp ~ -1900 km/s in H/3 and 
Ha lines alone. In subsequent spectra (day >21), we iden¬ 
tify a broad and structured absorption at Vexp ~ -1600 km/s 
and two additional components at ~ -550 and -750 km/s. The 
evolution of the observed P-Cygni absorptions of H/3 and Na 
I D2 is shown in Pig. H 

We have identifieoT via cross-correlation procedure be¬ 
tween observed narrow features and laboratory wavelength^ 
319 (over a total of 424 absorptions at day 7) low-excitation 
(Ei„ < 6 eV) transitions of singly-ionized heavy elements 
(Ba, Cr, Ee, Mn, Sc, Sr, Ti, V, and Y - see for example 
Williams et al.||2008)l, all of them characterized by an ex- 
panding velocity of Vgxp ~ - 550 km/s, see Pig|3 This pro¬ 
cedure resulted in some non-identifications andaegenerate- 
identifications, which we estimate to represent a small frac¬ 
tion (less than 20%) of the entire dataset. A complete analysis 
will be published elsewhere (Izzo et al. in preparation). 

Although from day 11 an additional expanding component 
is detected in the P-Cygni profiles of the Na I doublet, see 

^ http://lucagrb.altervista.org/research/lightcurve_NCen_I.pdf - courtesy 
AAVSO 

^ We referred to the Atomic Line List (v 2.05) maintained by Peter van 
Hoof, http://www.pa.uky.edu/~peter/newpage/ 
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Figure 1. The evolution of the radial velocity systems observed through the P-Cygni profiles in H/3 (left panel - centered at \np = 4861.32) and Na I D2 (right 
panel - Ajva/ = 5889.93) lines in the first three weeks of the V1369 Cen outburst. Expansion velocities of -550 and -1100 (±30 km/s) are marked with gray 
rectangles. Note the appearance of multiple expanding systems with increasing time. In particular, the absence in the spectra of day 4 and day 7 of the component 
at higher velocities {v^xp = —1350 km/s) in the Sodium P-Cygni profile, which is however observed in H/3. 



Figure 2. Spectra of V1369 Cen obtained on day 7, 13 and 21 between 4190 
and 4310 A with the identifications of some narrow absorptions lines. 


right panel in Fig. [u no further expanding components are 
identified for these heavy elements absorbing systems: the 
cross-correlation method provides a minor number of identi¬ 
fications considering a unique expanding velocity of ~ 1100 
km/s. If we consider the presence of both expanding com¬ 
ponents for these heavy elements, the number of acceptable 
identifications for the higher velocity components is a small 
fraction ~ 1% of the identified lower velocity systems. 

Among the many narrow transitions, we have identified 
many low-ionization neutral elements transitions (Fe, Ca, K), 
all belonging to the one expanding absorption system at the 
same velocity that observed in ionized heavy elements, with 
the only exception being Na I D lines, which show additional 
expanding components from day 11, see Fig. [T] In particular, 
we note the clear presence of the resonance transitions of ^Li 
I 6708, Ca I 4227, and K I 7699, see Fig. all of them with 
an expansion velocity of Vexp ~ -550 km/^ in the spectra of 
day 7 to day 18. 

2.1. Possible alternative identifications 

Here below we discuss several plausible alternatives to our 
^Li identification. 

1. We can exclude the possibil ity that this feature is due to a 
diffuse interstellar band (DIB - |Herbig|1995[[Bondar|2012] l, 
because i) no known DIBs are located at 6695.6 A; ii) all 
observed heavy elements narrow absorptions vanish after ~ 


20 days from the initial outburst when the nova continuum 
is still bright and DIBs should persist, being related to inter¬ 
stellar medium located in between the background exciting 
radiation and the observer; iii) we observe variations in the 
observed wavelengths of all these narrow absorptions, an ef¬ 
fect that DIBs do not show. 

2. An other possibility that must be considered is whether 
the 6695.6 A absorption line is due to metal lines excited by 
resonance absorption of UV radiation which is absorbed in the 
iron curtain phase and then rep rocessed at longer wavelengths 
(Johansson 1983[ Shore|20l2 1 . For example, each of the reso¬ 
nance transitions of Li, Ca, and K could be alternatively iden¬ 
tified with a transition of Fe or an s-process element at that 
wavelength, assuming that they are pumped by UV radiation 
(either continuum or lines). Particularly for the Li I 6707.8 
A resonance doublet, a possible identification could be Fe 
II 6707.54 A. In this case one should also expect absorp¬ 
tion features from the same lower energy level configuration 
(3d^CD)5p), and from the same term (4Fb-4G), which re¬ 
sults in additional 2 other possible transitions (Fe II 6769.27, 
6811.49 A) in the observed spectra. We have checked fo: 
their presence in the spectra of day 7 and day 1^ but we dii 
not find any absorption feature corresponding to Fe II 6769.27 
and 6811.49 A. This evidence disfavours an UV-pumped ori¬ 
gin for the absorption at 6695.6 A, even if considering the 
coupling between the transitions and the ejecta velocity field. 
Conversely, it supports our initial identification as expanding 
Li I 6708 A. 

3. We compare the red spectral region of Nova Cen with 
the same region observed in GW Or i, a T Tauri star charac- 
terized by the presence of Lithium (|Bonsack & Greenstein| 
|1960| l. GW Ori was observed with F'EROS and we have se- 
lected the observation of 26 November 2010. In Fig. |^we 
show the comparison of these two spectra, after correcting the 
nova spectrum for the expanding velocity of Vexp = -550 km/s. 
We clearly see the presence in the GW Ori of broad Ca 16718 
A and Li I 6708A giving more support to the identification 
of absorptions at 6695.6 A and 6705.3 A detected in V1369, 
with Li I 6708 and Cal 6718. 


^ we consider a A A = | A^j , — A/^f, | < 0.4 A. 

® these epochs corresponds to the first two EEROS observations, which 
show the largest number of narrow absorptions, due also to a greater resolu¬ 
tion, in our entire spectral database 
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Figure 3. The identification of Li I 6708 (left), Ca I 4227 (medium) and K I 7699 (right) features by direct compaiison with the Na I D2 5890 A P-Cygni 
absorption in the day 7 (upper panels) and day 13 (lower panels) spectra. All the features share same expansion velocity {Vgxp.Fi = —550 km/s, Vexp.Fi = —560 
km/s) as that of Sodium. 



Figure 4. The compaiison between the spectra of the VI369 Cen and GW 
Ori. It is clearly evident in both spectra the presence of ^Li I 6708 A absorp¬ 
tion line, as well Ca I 6718 A. The absorption in the nova spectrum around 
A = 6711.3 is identified as Cr II 6711.29 A. 


3. RESULTS AND DISCUSSIONS 

Under the assumption that this ab sorption is ^Li, w e have 
estimated its ejected mass following |Friedjung| (| 1 979|) . Since 
lithium, sodium and potassium are alkali metals with very 
similar Grotrian diagrams and with respective resonant tran¬ 
sitions, differing by only 0.25 eV for K-Li and Li-Na, we 
can assume that resonance doublets form under similar con¬ 
ditions. This assumption implies that the ratio of their optical 
thickness r, is related to their abun dance ratio mu ltiplied by 
their respective gf ratio. Following Spitzer|1998 (Eq. 3-48) 
we have for the case of Li/Na : 


Am(Li) _ / Wf,;6708 , ) gfNaDl Uu 

A^iNa) ~ 67082 ' 5890 ^) "" gfum% "" 7^’ 

where W\. the measured equivalent width at the transition 
wavelength A,, and m, the atomic mass of the corresponding 
element (in our cases un = 7, u^a = 23 and = 39). For 
the ^Li I doublet blend, we have considered the value of log 
gf = -0.174, whereas the single components has respectively 
lo g gf D1 = -0.00177, and log gf D2 = -0.3028 (jKramida 
et al. 2013 1 . We have obtained A„,(Li)/A,n(Na) = 3.2/100 
and Am{U)/Am(K) = 6.5/100 on day 7, and A,„(Li)/A^iNa) = 


2.4/100 and Am(T0/^m(^) = 5.2/100 on day 13. The lithium 
log overabundances are 4.8 and 3.9 with respect to s odium 
and potassium solar abundances ( [Fodders et al.|2009 1 . This 
result implies an overabundance of lithium in the nova ejecta 
of the order of 10 "^, which i s largely enough for explaining 
the Galactic ’Li enrichment ( Friedjun^|l979 1. The mass of 
sodium and potassium ejected in novae can be computed in 
terms of solar mass by using the results of different nova 
composition models, characterized by different WD masses , 
accretion rates and mixing degree pose & Hemanz||1998||. 
We have considered the results obtained for CO nova models, 
where the mass of ejected sodium ranges from 3.4 x 10 “^ to 
2.0 X 10“^ the mass of ejected hydrogen, while the ejected 
potassium varies from 5.1 and 7.2 x 10“®. The hydrogen 
ejected mass can be estimated both from the intensity of 
when the e jecta is completely optically thin, e.g. in late nebu¬ 
lar phase s ([Mason et al.|2005]l, a nd also through the observed 
t 2 value ( Della Valle et al.|2002[ ). Both methods converge to¬ 
ward the value of an ejected hydrogen mass of Mfj « 10 “"^ 
Mq. After combining the Na and K ejecta measurements with 
the respective lithium mass abundance ratio, we obtain that 
the mass of lithium ejected by VI369 Cen to be in the range 
Mu.i = 0.3 - 4.8 X10“'° M 0 . The ratios Li/Na and Li/K have 
been determined on day 7 an 13 when only a small amount 
of ’Be has decayed to ’Li. Hence, the measured amount of 
’Li may be only a lower limit. However, TNRs which pro¬ 
duce ’Be could s tart even years or mon ths (see Starrfield, II- 
iadis and Hix in [Bode & Evans[[2008| l before the nova was 
discovered. Therefore the most plausible scenario for VI369 
Cen is the one in which TNRs started weeks/months before 
the optical detection, and therefore the measured abundance 
of lithium can be considered a good approximation of the 
total amount of lithium actually produced by the nova (or a 
firm lower limit). The crucial quantity needed to compute the 
global Li yield of the galactic nova population is the Galactic 
nova rate, which is known within a factor two, = 20-34 
events/yr ( Della Valle & Livio|1994[jShafter|1997[ ). Wit 
e derive for the lithium mass ini 


With the 
injected 


yield obtained above, we 
in the Milky Way by nova systems Mu, tot ~ 2 - 45 M 0 /Gyr. 
However, it is well known that the ejecta of ’’slow” novae are 
more massive than ejecta of ’’fast” n ovae, by an order of mag - 
nitude, ^ 10“"' M 0 vs. ~ 10“^ M 0 (Della Valle et al. 2002 1 , 
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therefore fast novae, which can form ^ 30% of the nova pop - 
ulation of the Milky Way (jDella Valle & Duerbeck||1993|l, 
should contribute only marginally to the global ^Li yield. The 
above reported range of lithium mass decreases to Mu.tot <17 
M 0 /Gyr, for a rate of ’’slow” novae of 15-24 events/year, in 


good agreement with the theoretical predictions (Jose & Her- 
|nanz|1998). 

In Fig. I^it is shown the A(Li) vs. [Fe/H] observed re¬ 
lation compared with Galactic chemical evolution model re¬ 
sults. The Galactic chemical evolution model used is an up¬ 
dated version of t he chemical evolution model of IRomano 
et al.|(|1999| |2001[), th at is based on the two-infall model of 


Chiappini et al. (1997 1 , and that includes AGB stars (Karakas 
2010|l, super-AGE stars (Doherty et al.|2014ai. Galactic 


mic rays ( Lemoine et al.||1998| l and novae, as Li producers. 
In this model, the Galactic inner halo and thick disc form 
by accretion of gas of primordial chemical composition on 
a short timescale 1 Gyr). The gas is efficiently turned into 
stars as long as its density is above a critical threshold, below 
which the star formation stops. The thin disc forms out of a 
second episode of infall of gas of mainly extragalactic origin 
on longer timescales (7-8 Gyr in the solar neighborhood) and 
with lower star formation efficiency. It is worth emphasizing 
that the adopted mass assembly history is consistent with what 
is obtained for Milky Way-like galax ies in a full cosmologi¬ 
cal framework ( |Colavitti et al.|[20^] l. The black continuous 
line in Fig. [^is the best fit to the data obtained with a model 
with all Li sources, starting from a primordial Li abundance 
of AiLi) = 2.3 and by assuming that each "slow" nova (the 
adopted current slow novae rate is 17 events/yr) ejects on an 
average Mu = 2.55 x 10“'° Mq in agreement with the mea¬ 
surement of Mu = 0.3 - 4.8 X 10“'° Mq presented in this pa¬ 
per. The black dashed line shows the predictions of the same 
model when a high primordial Li abundance is adopted (see 
the Introduction). The red line is the best model with all Li 
factories but novae; it is clearly seen that novae are necessary 
to explain the late rise from the plateau value. The grey area 
indicates the uncertainties in the model predictions due to un¬ 
certainties in both the estimated Li yield from novae and the 
current slow nova rate; the upper (lower) boundary refers to 
the upper (lower) limit to the Li yield estimated in this pa¬ 
per and a maximum (minimum) current slow nova rate of 24 
(15) events/yr. The light green area similarly indicates the un¬ 
certainties in the model predictions when the m aximum and 
minimum Li yields from|Jose & Hernanz|(| 1998)1 are assumed. 
Though the two areas partly overlap, it is clearly seen that the 
theoretical nova yields tend to underproduce Li in the Galaxy, 
while the semi-empirical yields estimated in this paper give 
a better match with observed data points. Should the result 
presented here be confirmed by further observations of ^Li, 
Classical novae would stand as one of the major Li producers 
on a Galactic scale. 
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